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We argue that photon counts in a superconducting nanowire single -photon detector (SNSPD) are caused by 
the transition from a current-biased metastable superconducting state to the normal state. Such a transition is 
triggered by vortices crossing the thin and narrow superconducting strip from one edge to another due to the 
Lorentz force. Detector counts in SNSPDs may be caused by three processes: (a) a single incident photon 
with sufficient energy to break enough Cooper pairs to create a normal-state belt across the entire width of the 
strip (direct photon count), (b) thermally induced single-vortex crossing in the absence of photons (dark count), 
which at high-bias currents releases the energy sufficient to trigger the transition to the normal state in a belt 
across the whole width of the strip, and (c) a single incident photon of insufficient energy to create a normal-state 
belt but initiating a subsequent single- vortex crossing, which provides the rest of the energy needed to create 
the normal-state belt (vortex-assisted single-photon count). We derive the current dependence of the rate of 
vortex- assisted photon counts. The resulting photon count rate has a plateau at high currents close to the critical 
current and drops as a power-law with high exponent at lower currents. While the magnetic field perpendicular 
to the film plane does not affect the formation of hot spots by photons, it causes the rate of vortex crossings (with 
or without photons) to increase. We show that by applying a magnetic field one may characterize the energy 
barrier for vortex crossings and identify the origin of dark counts and vortex- assisted photon counts. 



PACS numbers: 74.78.-w, 85.25.Pb 

I. INTRODUCTION 

The superconducting nanowire single-photon detector 
(SNSPD) consists of a thin and long meandering supercon- 
ducting strip carrying a bias current / slightly below the criti- 
cal current I c , where the superconducting state is metastable. 
When a photon interacts with the strip it creates a hot spot 
in the film that drives a belt-like region across the width of 
the strip to the normal state. Consequently, a voltage pulse 
caused by the current redistribution between the supercon- 
ducting strip and a shunt resistor parallel to the strip is de- 
tected on nano-second time scales. After the normal belt of 
the strip cools down, the strip returns into the metastable su- 
perconducting state. Measuring voltage pulses, single pho- 
tons can be detected and counted. However, similar pulses are 
also recorded in the absence of photons, so-called dark counts, 
which introduce uncertainty in the counting of single photons. 

Recently, thermally activated single-vortex crossing was 
proposed as a possible mechanism for dark counts in 
SNSPDs. 1 As a consequence of the Lorentz force acting on 
a vortex crossing a thin and narrow current-biased strip the 
energy QqI /c is released, which for currents / > 0.6/ c suf- 
fices to create a normal belt across the entire width w of the 
strip (extending to a few correlation lengths £ along the strip). 
This process causes the transition from the current-carrying 
superconducting metastable state (S) of the strip at / > /* to 
the state with the normal (N) belt, and induces a current redis- 
tribution into the shunt accompanied by a measurable voltage 
pulse in the SNSPD. Here, I c is the critical current at which 
the energy barrier vanishes for vortex crossing. 

It was shown in Ref. [T] that in a thin, d < £ <C A, and 
narrow, w <C A, film the single-vortex crossing has the low- 
est energy barrier for creating dissipation, while phase slips 



across the entire width of strip and vortex-antivortex depair- 
ing are characterized by higher barriers. Here A = 2X 2 /d 
is the Pearl length with London penetration depth A and film 
thickness d. In fact, calculations accounting for the mecha- 
nism of vortex-antivortex binding and unbinding resulted in 
an energy barrier at least twice as large as for a single-vortex 
crossingPEI Since all these processes are thermally activated, 
that is, their rates depend exponentially on the barrier height, 
one can safely consider only the process with the lowest bar- 
rier. The vortex-crossing rate derived in Ref. [T] as a function 
of /, w and temperature T is in agreement with dark count 
rates measured by Bartolf et alP Recently, Hofherr et a/Palso 
ascribed the single- vortex crossing mechanism to the origin of 
dark counts in SNSPD. 

Here, we show that a weak magnetic field perpendicular 
to the superconducting strip results in a significant suppres- 
sion of the energy barrier for vortex crossing. This leads to 
an increase in the vortex-crossing rate. Hence, the vortex- 
based mechanism for dark counts can be tested experimen- 
tally. Also, we consider the case when the energy of a single 
photon is not large enough to create the normal belt across 
the entire strip width. Rather it leads to a belt-like supercon- 
ducting region with an elevated quasiparticle temperature. In 
this situation a photon-induced hot spot enhances strongly the 
probability of a subsequent vortex crossing. The combined ef- 
fect of hot spot and vortex crossing leads to the creation of a 
normal belt across the strip, i.e., to a single photon count. We 
call such photon counts "vortex-assisted". 

In recent SNSPD experiments, Hofherr et al 4 showed for a 
NbN strip, 97 nm wide, that photons of 500 nm wavelength 
give rise to a photon count rate that decreases at low cur- 
rents. Similar results were obtained by the NIST group^for 
a-W^Sii-^ based SNSPD. Our interpretation of these exper- 
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FIG. 1: (Color online) Cartoon of a superconducting strip with a 
vortex crossing without and with photon. Top: from left to right, 
illustration of the thermally excited vortex crossing and subsequent 
formation of a normal-state hot belt across the strip width resulting in 
a dark count. Bottom: an incident photon creates a hot spot (hot belt) 
across the superconducting strip, followed by the thermally induced 
vortex crossing, which turns superconducting hot belt into the normal 
state resulting in a vortex-assisted photon count. 



iments fits the situation described above, i.e., a single photon 
cannot turn an entire belt into the N state and additional en- 
ergy due to vortex crossing is needed for an S-N transition 
to happen in the belt. As the probability for vortex crossing 
depends strongly on the bias current and drops significantly 
at currents well below I c }^ the photon count rate becomes 
current dependent at low currents. 

Hence, our general picture is that a vortex crossing may 
trigger the S— >>N transition. A photon makes this process 
much more probable by creating a spot with suppressed super- 
conducting order parameter and thus with lower energy barrier 
for vortex crossing. A sketch of the strip and of the belt across 
are shown in Fig. [I] Like the bias current, an applied magnetic 
field enhances the rate of vortex crossings and causes an in- 
crease of the photon count rate. Therefore, the mechanism of 
vortex-assisted photon counts can be verified experimentally 
because the magnetic field does not affect the creation of hot 
spots, but increases the probability of vortex crossings. 



II. EFFECT OF THE MAGNETIC FIELD ON THE 
VORTEX-CROSSING RATE 

A. Basic equations 

In the following r = (x, y) is a position on the thin strip, 
< x < w, —L/2 < y < L/2, and w <C L. In the presence 
of a vortex at position r v , the superconducting order parame- 
ter, normalized to its uniform zero-current value, reads 

*(r,r v ) = /zexp{z[<p(r,r v ) - ny/t + ipo]}, (1) 

see, e.g., Refs. 7 8 Here fi 2 = 1 — n 2 describes the suppres- 
sion of the order parameter by the bias current. In Eq. ([I]), 
is a constant, i.e., its behavior in the vortex core is dis- 
regarded. Hence, we employ a variation of the London ap- 
proach in which, however, the Ginzburg-Landau suppression 
of the order parameter by transport supercurrents is taken into 
account. The parameter /i does not change much when the 
current varies from to I c , because 1 < l//i 2 < 1.157. The 



critical current I c is defined by a vanishing energy barrier, see, 
e.g., Ref. 1 , and is given by 
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e ~ 2.718. In this expression the parameter £ originates from 
the cutoff of the supercurrent at distances of the supercon- 
ducting correlation length in the London approach. Hence, 
the critical current, Eq. ([2]), determined in the London ap- 
proach has an uncertain numerical factor of order unity. Note, 
that the critical current depends on the film configuration, for 
more details see Ref. [9| It is easy to see that this current 
is of the same order as the depairing critical current density 
c<I>o/12\/37r 2 A 2 £ within Ginzburg-Landau theory, when mul- 
tiplied by the strip cross-section wd. Since is suppressed 
by the supercurrent 7, the Pearl length should be renormalized 
accordingly A = A/fi 2 . 

In narrow strips, w <C A, supercurrents are found by ne- 
glecting the induced magnetic field and the corresponding 
vector potential. 1 Then the sheet current density g for the vor- 
tex at position r v — (x V: 0) is g = curl Gz, 
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where X = ttx/w and Y = iry/w are dimensionless coordi- 
nates. The superconducting order parameter phase 
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Note that the characteristic length of the sheet current and 
phase variations for noninteracting vortices in both x and y 
directions is w. 1 In a strip of length L ^> w, a single vortex 
far away from the strip ends creates the phase difference 



<p(L/2)-<p(-L/2) = 2X v , 



(5) 



so that a vortex crossing from X v = to X v = tt changes the 
phase difference by 2tt. This is a generalization of phase slips 
in one-dimensional (ID) wires. The vortex crossing from one 
edge to another then causes the phase difference at the strip 
ends to vary in time. This induces a voltage change along the 
length L of the strip. The Josephson equation relates voltage 
and phase, which in our case relates the induced voltage with 
vortex velocity: 10 
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The energy of the vortex in zero applied magnetic field 
and at zero bias current is ($>o/2c)G(x — > x v ,0), whereas 
in the presence of the bias current I one needs to add the 
term <$>qIX v /7t due to the Lorentz force. In the applied field 
H = Hz, we also need to add the magnetic term —M V H, 
where M v is the magnetic moment of the vortex: 11 
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FIG. 2: (Color online) Position dependence, X v — 7TX V /w, of the 
potential energy for vortex entry in zero magnetic field for several 
bias currents. The potential is cut off for distances less than one 
coherence length £ from the left and right edges of strip, because of 
the finite vortex core. 



Finally, the potential energy U of a vortex in the presence of 
the bias current / and the external field H is given by: 



e^U(I, H,X V ) = In — smX v )--X v 
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Here eo = ^olo/nc = <£q/87t 2 A is the characteristic vortex 
energy in a thin film and Ho = <$>o/2w 2 . The potential energy 
in Eq. ^ corresponds to the vortex aligned with the applied 
field in agreement with Refs. 1111121 However, it differs by 
a factor of 2 in the logarithm with the potential employed by 
Gurevich and Vinokur. 6 Note that the equation ^ for the po- 
tential energy is valid only when the vortex core is far enough 
from the edges. 

In the following, we consider the situation of magnetic 
fields H < H c i(w), where H cl {w) is the field that corre- 
sponds to the vortex at X v = tt/2 (the middle of the strip). 
Thus the lower critical field for vortex entry is determined by 
£7(0, H cU tt/2) = and yields 



H cl (w) 
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In Figs.[2]through|4]we present the vortex energy as a function 
of position on the strip for different bias currents and magnetic 
fields. 



FIG. 3: (Color online) Position dependence of the potential energy 
for vortex entry for H — Ho for several bias currents. 
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FIG. 4: (Color online) Position dependence of the potential energy 
for vortex entry for H = 4Hq for several bias currents. 



B. Barrier for vortex crossing in magnetic field 

The energy barrier e&(J, H) is given by the maximum of 
U(x v ) in the interval £ < x v < w — £, where the expres- 
sion ([8} holds (For short we dropped the dependence of U 
on / and H and use the coordinates x v and X v interchange- 
ably). The positions of the extrema of U (x v ) are determined 
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FIG. 5 : (Color online) Graphic solutions for barrier heights for sev- 
eral currents % = I/I c o — 0, 0.1, 0.5 and fields h = H/H = 0, 1,4 
obtained by equating the LHS and RHS of Eq. ftp) . 
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FIG. 6: (Color online) The position X s of the maximum of the vortex 
potential U(x) (top) and the barrier height €b (bottom) as functions 
of applied perpendicular magnetic field for different bias currents. 
The thin lines with "+" symbols are the corresponding approximate 
barriers from Eq. (T2) , which are in excellent agreement with the 
exact solutions for larger currents. 



by dU(X v )/dX v = 0: 
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In Fig. [5] we plot separately the left- and right-hand sides 
(LHS and RHS) of this equation for different bias currents 
I and fields H to show qualitatively the solutions of Eq. ( [TO} . 
For H < ttHq 1 2 only one intersection exists at X v = X s cor- 
responding to the barrier energy e&. For fields H > ttHo/2 
and small currents one finds three intersections corresponding 
to a global maximum, local minimum, and local maximum of 
U (x v ), respectively, see Figs.|4]and[5] On the other hand, for 
large currents and large fields only one extremum (maximum) 
exists. At I = a local minimum exists in the middle of the 
strip at X v = tt/2 for fields H > 3Hq/2. However the vortex 
should overcome the barrier to reach this local minimum. The 
saddle point position for this barrier is shown in the top panel 
of Fig. [6j while its suppression with magnetic field is shown 
in the bottom panel. 

In the following, we discuss the case of I ^> Iq as Io is 
rather small in comparison with the currents on the order of I c 
employed usually in SNSPD experiments (I c ~ (w/£)Io ^> 
Jo). The position of the saddle point (global maximum) X s <C 
1 is given by 
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when X s > tt^/w. Thus the energy barrier 65(1, H), see also 
Fig.[6| for vortex entry is given by 



65 2wJo/(7reQ 

eo ~ n I + I0H/H0 ~ n /+ 



(12) 



where I c q is the zero-field critical current given by Eq. (2). We 
see that for large currents, one can account for the field effects 
by renormalizing the current: J + = I + IqH/Hq. Therefore, 
at high currents, the vortex-crossing rate R V (I, H, T) can be 
obtained from the zero-field rate R V (I, T) by replacing I — > 
1+ . This approximation holds for large currents not too close 
to J c o to fulfill the condition X s > tt^/w. From Eq. (12) it 



follows that for sufficiently small magnetic fields, H <C H* 
the field dependence of the critical current is linear in H, 



Ic(H) = J* 1 



H 
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C. Vortex-crossing rate 

In Refs. 11161 the vortex-crossing rate was obtained by em- 
ploying the known stationary solution for a potential U(x v ) 
with periodic boundary conditions (i.e., for the potential 
U (x v ) created by periodic extension to the whole axis — 00 < 
x v < +00]. Here, we derive the vortex-crossing rate for the 
realistic potential U(x v ), which is not periodic. In order to 
correctly treat vortex-assisted photon counts, we need a time- 
dependent solution. We will use the rate derived in the quasi- 
stationary approach to find vortex crossing rate via relaxing 
hot spot and we will present condition when such an approach 
is valid. Besides, we show that the dissipation in the shunt 
resistor is comparable to the Bardeen- Stephen dissipation in 
the vortex core and cannot be neglected considering the dissi- 
pative vortex dynamics in SNSPDs. 
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1. Langevin equation 

The Bardeen-Stephen dissipation power is rj v Xy, where 
f] v = <3>o/(27r£ 2 c 2 i?n) is the Bardeen-Stephen drag coeffi- 
cient. Here R\j is the film sheet resistance taken slightly above 
T c . The power dissipated in the shunt with resistance R s , 
due to the voltage pulse of a crossing vortex, is V(t) 2 /R s = 
(<$>ox v /cw) 2 /R s = f] s xl, see Eq. ^ for voltage-phase rela- 
tionship. Therefore the equation of vortex motion in the pres- 
ence of shunt and thermal noise is 



dU(x v ) 
dx v 



F L (t), 



(14) 



where Fl is the Langevin force with the effective drag coeffi- 
cient 



Veff = r) v +r) 8 
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For the parameters w = 100 nm, £ = 4 nm, d = 6 nm, 
i?n = 400 £1 and R s = 9 ft the shunt contribution to the 
viscosity is « 0.45/^ and thus cannot be neglected. 



2. Fokker-Planck equation 

To find the vortex-crossing rate for diffusive vortex motion, 
we use the Fokker-Planck equation and Kramers' procedure 
to calculate the escape rate for particle diffusion over a bar- 
rier as described in Ref. 13 As mentioned previously, in the 
London approximation, one can treat the vortex as a particle 
described by the potential energy in Eq. ^ for coordinates 
xo < x v < w — xq. At the left edge, i.e., at distances smaller 
than xq ~ £ an excitation precedes the crossing of the vortex 
(so-called "pre- vortex"). In the presence of the bias current 
with its self-induced magnetic field lines at the edges, the pre- 
vortex has a normal core, but no circular supercurrents en- 
closing it. Therefore, within the London approximation the 
boundary condition on the left side may be recast in the form 
of a trapped particle (pre- vortex) in a box potential sketched 
in Fig. [7] with the effective potential 



U(x v ) 



oo 




U(x v ) 



x v — 0, 

< Xy < X0, 

X v > Xq. 



(16) 



The infinite repulsive potential U at x v =0 prevents the 
pre-vortex from escaping to the left, because pre-vortex 
or vortex (particle) exist only inside the superconducting 
strip. The probability density of such particle is W(x v ) w 
x^ 1 exp[— U(x v )/T] « x^ 1 between < x v < xq, where 
we take into account the uncertainty of the vortex coordinate 
of the order xq. At x v > xo the vortex is fully developed in- 
side the strip and its diffusion over the barrier between points 
xq and w — xq is described by the potential U (x v ). 

Finally, we solve the Smoluchowski (Fokker-Planck) equa- 
tion for a weakly time-dependent potential U(x v ,t) to ac- 
count for changes of the superconducting order parameter 



with time due to the hot spot relaxation: 

dtW(x,t) = -d x S(x,t), (17) 
S(x,t) = -De~ (j ^ t ^ T d x \e V ^' T W{x,t)\ ,(18) 

where D = T/r] e ^ is the diffusion coefficient, W(x, t) dx is 
the probability to find the vortex in the interval (x, x + dx) at 
the time t, while S(x, t) is the probability current density. We 
consider bias currents I < I c , such that e^/T is large and the 
crossing rate is small. Then S and d t W are small, while the 
probability current density S(t) = S(x,t) is approximately 
independent of x, because d x S w 0. At the left edge, the 
probability density of a vortex is W(xo,0) = 1/xq. At the 
right edge, W(w — xq, t) is exponentially small as the vortex 
escapes and W(w — xo, t) = can be taken there. Integrating 
Eq. ( [T8] ) between x v = xq and x v = w, we obtain for the 
current probability 



rW 

(D/x )W(x ,t) = S(t) / dx< 

J Xn 



V(x,t)/T 



(19) 



Our solution for S(t)~ r coincides with the one found by 
Kramers 's escape rate for the double- well potential, see 
(5.109) in Ref. [H if we set the potential f(x) = in the 
interval (xi, X2) and take p = 1 in (5.109). Also neglecting 
the exponentially small probability for vortex crossing from 
the opposite side, that is, from x v = w to x v = 0, we find the 
vortex-crossing rate 



R v (t) = S(t) 
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X 



dxe u ^' T 



(20) 



Note that by symmetry the rate of anti vortex crossing from 
x v = w to x v = is the same as for vortex from x v = to 



w. The main contribution to the integral in (20) comes 
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FIG. 7: (Color online) The effective potential U(x v ) in the presence 
of bias current for a vortex in the particle approximation (red circle). 
The left- side boundary condition of the Fokker-Planck equation de- 
scribes the existence of a "pre-vortex" within distance xo ~ £ of the 
strip's left edge. 
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from regions near the maximum of U (x), where it coincides 
with U (x) for currents not too close to the critical one. Under 
this condition we can make the approximation 



R- X D « x / dxe u{x ^ /T 
Jo 



(21) 



Further, we can replace the prefactor x with £ to obtain 
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Z(tt), 



(22) 



where p = zaT/2o and v = co/T and eo(t) and Jq(£) are 
assumed weakly time-dependent in the potential of Eq. ([5]). 
The integral Z(7r) h as an exact solution with an asymptotic 
expansion for v ^> lJHIU 



7rexp(-7rp/2)r(;y + 1) 
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with the Gamma function T(z). Finally, for / close to I c q 
we find the total vortex and antivortex-crossing rate R tot = 
R v + R av — 2R V at H = in a belt of size w x w: 



47rTc 2 i? e ff 



J 



/co(*) 

27 



y(t) + l 



(24) 



xe 3 



ii 



eff 



(7re^ c0 (t)/2 W /) 2 

^ + 2^/w) 2 ifc ' 



(25) 
(26) 



The solution for the parameters changing with time on the 
characteristic time scale r the adiabatic approach is valid if 
the thermalization time is much shorter than r, i.e., if 
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This condition is satisfied if r <C 0.1 ps for R e fi = 400 ft and 
A = 60 /im. 

For currents / ~ J c o the exponential factor on the right- 
hand-side of the attempt frequency Q is of order ~ 1/e for 
materials parameters of interest (w/£ ~ 25 and ^ ~ 110). 
Further assuming i? e ff = 400 T = 5 K, and / = I c q, we 
estimate ^ ~ 4.4 x 10 10 s _1 . However, one should not put too 
much weight on the estimate of the 'pre-exponential' factor ft 
compared to (J// c o) zy+1 , because the approach "vortex as a 
particle" was used and such an approach is only approximate 
near the strip's edges, resulting in possible corrections of order 
unity. 

It is worthwhile to point out that our expression for the rate 
R v is smaller by a large factor T{y - l)w/£ = 3.5 x 10 36 



compared to the voltage-current result derived in Ref . 6 There 
are three reasons for this significant discrepancy: (i) the fac- 
tor 2 U originates from the logarithmic difference in the vortex 
potential U in Eq. ([8]), (ii) the factor (y — 1) is due to the 
difference in boundary conditions (periodic vs. non-periodic) 
for a finite-size geometry, and (iii) the factor w/£ stems from 
the fact that vortices are uncorrelated (non-interacting) at dis- 
tances larger than w along the y-axis and not £, as was as- 
sumed in Ref. [6| 

Another noteworthy result of this work is that a shunt 
suppresses the rate because the effective vortex viscosity in- 
creases due to the additional dissipation in the combined sys- 
tem of superconducting strip and normal shunt. Finally, we 
note that at large v the condition of small probability current 
density S(t) is fulfilled for practically all currents of inter- 
est. Thus justifying our original assumptions for solving the 
Smoluchowski equation. 



3. Field-dependent rate 

To evaluate the magnetic field effect on the rates, we replace 
/ with J + of Eq. ( [T2] ) as discussed above. Then, one can see 
that the vortex-crossing rate increases with field according to 



R v (I : H)=R v (I + )^R v (I,0) 1 



I c qH 
IH* 



• (28) 



For I c oH <C IH* and in the large limit for the exponent, 
v ^> 1, we obtain for the vortex-crossing rate 



R v (I,H)^R v (I,0)e H / H \ H x = 



H*I 



(i/ + 1) J, 



(29) 



cO 



Next, we consider the crossing rate for antivortices entering 
from the right edge x = w and moving toward x = 0. The 
zero-field potential barrier for antivortex is a reflection of that 
for a vortex in the strip middle, i.e., the zero-field rates for 
vortices entering at x = and for antivortices starting from 
x = w are the same. However, the field contribution —M a H 
to the antivortex energy has opposite sign to that of a vortex. 
It is straightforward to obtain the corresponding barrier 
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cO 
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(30) 



We then obtain the rate for antivortices after similar calcula- 
tions 



R av (I,H)^R v (I,0)e- H / H K 



(31) 



The total crossing rate is the sum of vortex and antivortex 
crossings 

R tot (I, H) = 2R V (1, 0) coshiH/H^ . (32) 



D. Characteristic magnetic fields 

Next, we estimate all characteristic fields for typical film 
and materials parameters of the SNSPD. We consider a film 
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with width w = 100 nm, thickness d = 6 nm, coherence 
length £ = 4 nm, and renormalized Pearl length A = 60 
/im at temperature T = 5.5 K with the parameter u w 110, 
see Ref. [T] In this case, the characteristic magnetic field for 
describing the potential change is 77 = 0.1 T, see Eq. ([5]). 
The field for vortex entry is H c \ = 0.36 T, see Eq. ([9]). The 
field scale for the critical current I C (H) is H* w 0.6 T, see 
Eq. ([131, and for the vortex-crossing rate is Hi ^4.1 mT at 



$o/(47r 2 A(i£ 2 ). We attain for the critical threshold current 



I = 0.757 c0 , see Eq. (29). 



III. TRANSITION S^N CAUSED BY VORTEX CROSSING 

It was argued in Ref. 1 that a single vortex, crossing at 
high currents, can destroy the superconductivity inside the 
belt, which covers the entire width of strip. If so, the vor- 
tex crossing results in the switching of the current from su- 
perconducting strip to shunt in SNSPD and the destruction of 
superconductivity in the whole strip, when the shunt is absent 
or insufficient. Here, we show that it occurs in strips with the 
width smaller than some critical w c . We will estimate w c and 
the current 7* above which it happens. 

A crossing vortex releases the energy ^qI/c as a cloud 
of quasiparticles is squeezed out of the moving vortex core 
as described by Larkin and OvchinnikovP^ These quasiparti- 
cles diffuse from the vortex path during crossing time To « 
w 2 $o I '(27r£ 2 cR\jI) of the order of several pico-seconds for 
w = 100 nm and bias currents of the order 7 C . This estimate 
was obtained in Ref. 1 assuming that the supercurrent is not 
affected by vortex crossings and thus the vortex moves under 
the effect of almost constant Lorentz force. During the time 
To of the order of several pico-seconds, the energy transfer to 
phonons and substrate is negligible.^ Thus, when the vortex 
exits, the quasiparticle cloud has approximately the form of 
a belt of width £(I) « (Dr ) 1/2 , where D is the diffusion 
coefficient. We estimate the width as 



1(1) = t 



( 2enAwD I c 



1/2 



4.29£ 



1/2 



(33) 



where R\j — 400 ^ is the film sheet resistance, D « 0.46 
cm 2 /s is the diffusion coefficient, and A = 60 /im, £ = 4 nm, 
w = 100 nm. We estimated the diffusion coefficient from 
D = 7r 2 kg/(3q 2 jdRc\), where q is the charge of the elec- 
tron, ks is Boltzmann's constant, and 7 = 220 J/m 3 K 2 is 
the Sommerfeld coefficient. The energy released by a vortex 
crossing, <&qI /c, suffices to destroy superconductivity in the 
area of the quasiparticle cloud (belt) at currents I > 7* , where 
7* in zero-field is defined by the energy balance 



7* _ el(P) _ ( e 3 AwD 
T c ~ 8tt£ ~ I 32ttc 2 R b £ 3 



1/3 



0.60. 



(35) 



To have 7* < 7 C one needs to satisfy the condition for the crit- 
ical width of the strip with materials parameters used above 



w c < 1.745 • 10 b r/A, 



(36) 



otherwise a single- vortex crossing does not result in the for- 
mation of normal belt across the entire strip width and the 
superconducting state is stable with respect to single-vortex 
crossing. For films studied so far, we have A/£ ~ 1.5 • 10 4 



and thus the condition Eq. ( 36 ) is fulfilled for strips of widths 
up to w c « 116£. In this case, the vortex crossing destroys 
the metastable state at currents 7 > 7* without shunt and pro- 
duces a voltage pulse in the combined system with shunt. Our 
estimate for 7* > 0.67 c is rather crude. It neglects the possi- 
ble creation of a normal region behind the moving vortex and 
the corresponding change of the supercurrent density in front 
of it. 



IV. VORTEX-ASSISTED PHOTON COUNTS 

Absorption of a single photon in a thin and narrow super- 
conducting strip results in a cloud of quasiparticles with en- 
ergies above the superconducting gap. First, this cloud of ex- 
citations grows in number of quasi partic les and in size due to 
avalanche processes and diffusionpE^ but in time the cloud 
diminishes and finally vanishes due to relaxation processes. 
We assume that in sufficiently narrow strips, w < 100 nm, (a) 
the cloud covers the entire width w of strip when it reaches 
its maximum size, (b) the quasiparticle density in the cloud 
(hot belt) is close to uniform, and (c) quasiparticles suppress 
the superconducting order parameter inside the hot belt, but 
their density is not sufficient to convert the hot belt to the 
normal state. Thus the superconducting condensation en- 
ergy density in the hot belt, 77 2 ^/87r, satisfies the inequali- 
ties < 77 c \/8tt < 77 c 2 /8tt. Further, we assume that quasi- 
particles in the hot belt thermalize exponentially with relax- 
ation time r of the order of 40 ps for film thickness consid- 
ered, see Refs. 1191201 It follows that the vortex-crossing rate 
via hot belt (denoted by subscript h), R v h(I), is enhanced 
in comparison with that for dark counts, because the param- 
eter Vh = coh/T ~ 77 2 /l /(87rT) is reduced in comparison 
with vq. In the framework of assumptions (a) and (b), the 



crossing rate inside the hot belt is given by Eq. (32), how- 
ever, with renormalized parameters eo^ < eo, diminished 
exponent Vh = e^jT < and reduced critical current 

Ich = Ico(vh/vo) < IcO- 



$0 



7* 



Hi 
8tt 



£(I*)wd. 



(34) 



Our estimate for the hot area results in a lower bound for 
the threshold current 7*, where the thermodynamic critical 
field of the superconducting condensate is defined by 77 2 = 



A. Current dependence of assisted counts 

There are three superconducting current-biased regions to 
be considered, which are separated by three characteristic cur- 
rents I I < I c h < ho, see Fig. [8] For currents below the hot 
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NO COUNTS CURRENT DEPENDENT RATE^ PLATEAU 



ch 



l c0 



FIG. 8: (Color online) Three characteristic current regimes and three 
characteristic behaviors for photon count rates. 



threshold, 1%, counts are absent, because vortex crossings do 
not result in the formation of a normal-state belt. 1% is de- 
termined through the energy balance for destroying the super- 
conducting condensate in the hot belt similar to that described 



in Eq. ( 34 ) for dark counts, however, with the thermodynamic 
field H ch ~ H c ^/v h /vQ < H c : 



(37) 



and hence 1% 



0.6I c h. Between 1% and I c h, the vortex- 
assisted photon counts are strongly current dependent with a 
precipitous drop below I c h- At high currents, I c h < I < 
I c o, photon count rates saturate at a plateau. In the region 
1% < I < I* ~ 0.6/ c o, vortex crossings via the hot belt 
destroy the superconducting state inside the belt during the 
relaxation time r of the hot spot. During this time, the belt is 
in the normal state, which causes the "vortex-assisted" photon 
count. 

If I c h < I < ho, the rate of single-photon counts, R pC9 is 
the same as the rate for hot spot formation Rh, because the 
barrier for vortex crossing is now absent. For hot counts Rh, 
we have Rh = r]hR p , where R p is the photon rate and r]h 
is the quantum efficiency of hot spot formation caused by a 
single photonPM 

Next, we consider the region Ifc < I < I c h, where the 
photon count rate depends on the bias current. This depen- 
dence has been observed in NbN strips with d = 5 nm and 
w = 97 nm for photons with 500 nm waveleng as well as 
in a-W^Sii-^ SNSDPs (d = 4.5 nm, w = 150 nm) for a wide 
range of photons of 672-1 850 nm. 5 In this current regime, the 
photon count rate is R pc = RhP^R^ 1 ), where V(t) is the 
probability of vortex crossing before time t; t is counted from 
the moment of the hot spot formation. For R pc , we take V(t) 
at t = R^ 1 , because a given hot spot does not contribute to 
photon counts if the vortex crossing does not occur before for- 
mation of another hot spot by the next incident photon. The 
probability V(t) satisfies the rate equation 



dV 

~dt 

which yields 

P(t) = 1 - 



= R vh (I,t)(l-P), 



exp 



■ / dt'R vh (I,t' 
Jo 



(38) 



(39) 



Due to the thermal relaxation of the hot spot the vortex- 
crossing rate via the hot spot changes with time, because 
I c h(t) and Vh(t) diminish with t. We will account only for 
time dependence of v as the crossing rate is much more sen- 
sitive to v than to I c h. For v(t) we have 




0.16 




W (0) 



FIG. 9: (Color online) The vortex-assisted photon count rate R pc vs. 
bias current given in Eq. ( |45] >. Top panel: Rates for several values 
of hot belt lifetime r at T = 5 K. Bottom panel: Rates for several 
temperatures at r = 1 ps. Other parameters for NbN are vk — 40 
and z/ = 110 at T=5 K, T c = 15 K, R p = R h = 10 6 s" 1 , R= 
400 Q, w — 100 nm, and £o = 4 nm. Inset: Current dependence 
of photon count rates in a-W x Sii_ x for different bath temperatures 
from Ref. |5] Note, that film temperature is generally higher than 
indicated due to heat released by photons and corresponding vortex 
crossings. 



where is is defined in the absence of the hot spot, that is, it is 
the same as for dark counts. Defining the integral in V(R'^ 1 ) 
in Eq. ([39]) as 



we obtain 



dtR vh {I,t) 



o 

dt 



dv + 1 — v ' 



. dt f v 
dv—R vh {I,v), 
dv 

(41) 



= vo + 1 + {yh ~ ^o) exp(-t/r), 



(40) 



V"h Jv h +1 v + 1 - V \IchJ 

J— I (42) 

where v\ = v[R h x ) is defined by Eq. ( |40| ). The integrand 
increases rapidly as v changes from Vh to a larger value v\ w 
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Vq + 1 when R^t <C 1. Hence, the main contribution comes 
from the region near Vh- We define v = 1/^ + 1+ a and expand 
the integrand in a. Integrating over a we find 



TR v (I,v h ) 1-(W/) 



Vq - Vh 



Wch/I) 



where 



R v (I,v h )^n(I,T, v h ) 



Lch 



Finally, we attain the vortex-assisted photon count rate 

R P c(I, v h ) =Rh[l- exp(-K(I, v h )}. 



(43) 



(44) 



(45) 



The measured count rate includes dark counts in the remaining 
"cold" part of the strip, which must be added to the right-hand- 
side of Eq. (45). Roughly, if R pc ^> R V (I, vq)L/w one can 



neglect the contribution of dark counts from the "cold" section 
of the strip. 

At high currents / > I c h, we obtain R pc « R p , that is, 
the photon count rate has a plateau as a function of /. Here, 
experimental data on photon rates provide information on the 
rate of hot spot formation by incident photons, R p . For lower 
currents in the interval I c h > I > 1%, the rate of photon 
counts shows a power-law current dependence, 

Rpc(I) * R h K « R h rR v (I, v h )/(v - v h ) (x I Vh+1 . (46) 

The low-current part of the curve In R pc vs In / allows us to 
extract Vh characterizing the suppression of superconductivity 
in the hot spot. A sharp crossover from high to low current 
behavior occurs at current values close to the critical current 
I c h due to high values of Vh- This crossover depends weakly 
on r, as one sees from the top panel in Fig. [9] Furthermore, 
at bias currents below 1% photons do not lead to formation of 
normal belts and thus are not counted. 

To prove that vortex crossings are involved in photon counts 
at / > II one may use the magnetic field which enhances the 
vortex-crossing rate along with photon counts at low currents. 
For H <C H*, and accounting for both vortex and antivortex 
crossings in the hot belt, this enhancement is described by the 
relation 

R pc (I, H) =R h (l- exp[-2ft( J, v h ) cosh(H / H lh )}] 

(47) 

with field scale H lh = #i(z/ + + 1). The field ef- 

fect is most pronounced at low currents when 11(1, Vh) <^ 1. 
Since the magnetic field affects vortex crossings, but not the 
formation of hot spots by photons, the dependence of vortex- 
assisted counts R pc on H will unambiguously confirm the in- 
volvement of vortices in photon counting. 

For wider strips, when the hot spots of excited quasiparti- 
cles do not cover all of the strip width, the hot spot diminishes 
locally the effective width of the superconducting strip and 
roughly one can replace w ^ w = w — D in Eqns. ( [24] ) 
and ([26]), where D is the hot spot diameter. Thus we can 
use Eq. ([47]) with the renormalized width w and the same vq 



as for dark counts. It is worth to remember that the vortex 
crossing rate depends strongly on w, R v (I,vq) oc w~ v °~ l , 
with v w 110. Hence even a small decrease in w enhances 
dramatically the vortex-crossing rate due to the formation of 
hot spots. Clearly, for a more rigorous treatment one needs 
to consider vortex crossings for a nonuniform superconduct- 
ing order parameter and nonuniform currents. However, our 
main result will stay: the magnetic field renormalizes the bias 
current as described above and results in the enhancement of 
photon count rates. 



B. Comparison with experiments 

In the framework of our model, the data presented in Fig. 4 
of Ref. 4 for nano wires with w = 97 nm yield an estimated 
Vh ~ 40 at 6 K. On the other hand, for strips of similar width 
and thickness, vq for the dark count rate at 5.5 K was esti- 
mated as vo « 110, see Ref. Q] where the data reported in 
Ref. [3] were used. Hence, the barrier suppression by the field 
is strong; still, more measurements to extract vq and Vh are 
needed, especially for narrow strips for which our model of 
uniform hot spots across the whole strip width is applicable. 

The part of the curve R pc (I) at intermediate currents pro- 
vides information on the relaxation time r. The dependence of 
Rpc(I) for Vh = 40 and vo = 110 at 5 K for several values of 
r and temperatures is presented in Fig. [9] We model the tem- 
perature dependent parameters I c h(T) = l c h(0)s^, £(T) = 
Zo/st, and Vi(T) = ^(0)4/T with e T = (1 - T/T c )^ 2 . 
Recently, similar behavior has been reported by the NIST 
groupS for the current dependence of the photon count rate 
in a-W^Sii-^ strip with T c ~ 3.8 K. The data for R pc /R p 
are shown in the inset of the lower panel of Fig. [9] 

Both calculated and observed current dependences of pho- 
ton count rates have a plateau at high currents and a sharp drop 
below the current associated with I c h- Note that the supercon- 
ducting condensation energy in the hot belt is proportional to 
^olch/c. Therefore one can anticipate that the condensation 
energy in the belt drops with increased photon energy. The 
wavelength dependence of the count rate in Fig. 2(b) in Ref. [5] 
is in agreement with such an expectation. 



V. DISCUSSION AND CONCLUSIONS 

In a somewhat oversimplified picture presented in the litera- 
ture, the hot spot formation results directly in a voltage pulse. 
The argument goes like that: The hot spot becomes normal 
and the transport supercurrent should redistribute in the re- 
maining superconducting regions. There, it exceeds the criti- 
cal value and transforms this area to the normal state as well. 
Because of the large resistance of the normal belt created, it 
induces the current redistribution to the shunt resistor and a 
voltage pulse is detected. To quantify this scenario, we pre- 
sented in this work a vortex-based phenomenological model 
for the transition from the current-carrying metastable super- 
conducting state to the current-carrying normal state. Two 
important points in our scenario are the metastability of the 
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current-carrying superconducting state with respect to vor- 
tex crossings and the energy balance of the transition process 
driven by such a crossing. 

One-dimensional wires (nanowires) with w < £ and two- 
dimensional wires with w ^> £, considered in this work, have 
many features in common. As we have already shown in 
Ref.[TJ phase slips in ID wires 7 8 play the same role as vortex 
crossings in 2D wires at low temperatures considered above. 
Each phase slip in nanowires releases the energy &oI/c. The 
critical current in nanowires, I c q = c$o£/(12v / 3 / 7r 2 A 2 ), is 
only slightly larger than for the strip. For currents / = I* > 
0.3/ c o, the energy released by phase slips is sufficient to de- 
stroy superconductivity in a segment of the nanowire of cross- 
section 7r£ 2 with length ~ 1.2£. Eventually, the whole wire 
becomes normal due to Joule heating if the shunt is absent. 
Therefore, in this case there is no dc voltage to be considered 
in the superconducting state. Similarly, at currents well below 
« 0.3/ c o the probability of phase slips is negligible and no 
dc voltage may be observed. On the other hand, in the pres- 
ence of a shunt, voltage pulses will be detected in full analogy 
with dark counts in shunted strips. Note, that the shunt resis- 
tance is even more important in phase slip rates than for strips, 
because the small factor (, 2 /w 2 in Eq. (26) of our model be- 
comes of order unity in nanowires and strongly renormalizes 
the viscosity. 

Interestingly, the proposed SNSPD mechanism is similar to 
that of the bubble chamber used to detect charged elementary 
particles in particle physics PD Inside the chamber, a super- 
heated liquid is prepared in a metastable state. On entry of a 
charged particle, the ionization track is formed around which 
the liquid vaporizes by forming mesoscopic bubbles. The bub- 
ble density around the track is proportional to the particle en- 
ergy loss. In SNSPD too, the transition from the metastable to 
stable states is used for the detection of a single photon. 

We summarize the calculations within the phenomenologi- 
cal model of vortex-assisted photon counts in narrow and thin 
SNSPDs by listing the main results: 

(a) We estimated the critical width w c of the superconduct- 
ing strip below which single-vortex crossings create a normal 
belt across the entire strip width resulting in dark counts in the 
SNSPD and transition to the normal state in the absence of a 
shunt. 

(b) We have shown that the magnetic field in combination 
with a bias current enhances the vortex crossing rate. At cur- 
rents / ^> Iq, close to I c , the effect of the magnetic field 
is equivalent to an increase of the bias current by a factor 



(1 + H/H*), where H* « 0.6T in 6nm thick and lOOnm 
wide films. The corresponding characteristic field affects the 
vortex-crossing rate (dark counts), because the effective cur- 
rent is much smaller, roughly by a factor 1/vo <^ 1. For 
films of similar geometry, fields H <C H* , and for currents 
/ = 0.75 ic(0) the vortex-crossing rate is renormalized by a 
factor cosh(iJ/.ffi), where E x ~ H*/v ~ 4mT. 

(c) We have derived the photon count rates that require a 
vortex crossing to transform the photon-created hot spot to a 
normal-state belt across the strip's width. A single photon in- 
cident diminishes the critical current in the hot spot eliminat- 
ing the barrier for vortex crossing at currents / > I c h (within 
the plateau region of the count rate). This results in a photon 
count when the photon energy is too small to create a normal- 
state belt. At currents below I c h a reduced barrier still exists 
and the photon count rate depends strongly on the bias current. 

(d) We also argued, based on available data, 4 5 that vortex 
crossings are the origin of photon counts in the SNSPD. This 
vortex-assisted mechanism may be verified by the application 
of magnetic fields, which enhances the vortex crossing rates, 
but does not affect the creation of hot spots by photons. We 
predict that already weak magnetic fields of order Hi should 
cause the enhancement of photon counts, in particular at low 
bias currents. 

Finally, our calculations show that the vortex-assisted pho- 
ton count rate is characterized by a plateau at high currents 
with a tail toward low currents. This feature suggests that an 
optimum current exists, which minimizes the effect of dark 
counts, while only weakly diminishing photon counts. We 
propose to operate SNSPD devices above I c h, the low-current 
edge of the plateau, where photon-induced counts are still ef- 
fective, whereas the dark count rate is strongly suppressed. 
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